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Abstract 



We present the program summary and long write-up for Sphinx tt, 
a new Monte Carlo simulation for transversely polarized nucleon- 
nucleon scattering. The program is based on Pythia/ Jetset. 



PROGRAM SUMMARY 



Title of program: SPHINX TT 
Catalogue number: 

Program available from: Oliver Martin, 

|http : / /th.pliysik.uni — trankturt.de/ ~ martin/sphinx.html 



Computer: IBM RS6000, PentiumPro 200 based PC; others with a FORTRAN 77 
compiler 

Operating system: tested under AIX and LINUX but does not depend on the 
particular operating system 

Programming language used: FORTRAN 77 

High speed storage required: 5 Mbytes 

Number of bits per byte: 8 

Number of lines in combined test deck: 19980 

CPC subprograms used: Jetset 7.3 

Keywords: transversely polarized nucleon-nucleon scattering, transversity 
parton distribution, high energy physics, Monte Carlo simulation 

Nature of physical problem: This program can be used to simulate trans- 
versely polarized nucleon-nucleon collisions at high energies. Spins of collid- 
ing particles are taken into account. The program allows to calculate cross 
sections for various processes. 

Method of solution: The existing Monte Carlo program Pythia 5.6 has 
been modified to incorporate spin effects. The program incorporates nearly 
all features of Pythia. 

Restriction on the complexity of the problem: The spins of the colliding 
hadrons must be transverse with respect to the collision axis. (A program 
for collisions of longitudinally polarized hadrons (Sphinx) is also available.) 
Furthermore, the spins of the two hadrons must be either parallel or antipar- 
allel with respect to each other. 

Typical running time: ~ 0.01 sec CPU per event on PentiumPro 200 based 
PC, depends strongly on kinematical cuts 

References: For general information see the Pythia manual [Q]. Specific de- 
tails on Sphinx TT (and Sphinx) can be obtained from: Oliver Martin, In- 
stitut fiir Theoretische Physik, Goethe Universitat, Robert-Mayer-Str. 8-10, 
60054 Frankfurt am Main, Germany, e-mail: martinOth . physik . uni-f rankf urt . de 



Long Write-Up 



1 Introduction 

During recent years the topic of transverse spin effects has become increas- 
ingly popular. Especially the so called transversity parton distribution 6q{x, Q"^) 
which is comparable to Ag(x, Q^) in the case of a longitudinally polarized nu- 
cleon has attracted much attention (for a concise introduction see [H). Due 
to experimental limitations, Q^) could not be measured so far. Only 
high energy collions of two transversely polarized nucleons or semi-inclusive 
deep inelastic electron-nucleon scattering with transverse polarization allow 
a determination of 5q{x, Q"^). 

However, in 1992 the RHIC Spin Collaboration was constituted 0. One of 
its major projects is the measurement of the transversity parton distribution 
by analysing collisions of transversely polarized protons at RHIC, which is 
now under construction at BNL. RHIC will be able to scatter protons at 
cm energies of up to y/s = 500 GeV with high transverse or longitudinal 
polarization of up to 70% . A first measurement may already be possible 
in the year 2000 or shortly after. To plan for such experiments requires 
extensive Monte Carlo simulations. Sphinx TT was developped to perform 
such studies. 

Developing a state-of-the-art Monte Carlo simulation takes several years. 
Therefore, we used Sphinx |^ which itself is an extension of the popular 
Pythia 5.6/Jetset 7.3 |T| as a basis for our program. Sphinx is able to 
simulate collisions of two longitudinally polarized nucleons in the framework 
of the QCD improved parton model and describes nucleons on twist-2 level 
[0]. On this level no new physics is involved in collisions of a transversely 
polarized nucleon with a longitudinally polarized nucleon ^ . For this reason 
we decided only to cover the scattering of two transversely polarized protons 
with our program, which we consequently called Sphinx tt. Sphinx is an 
acronym for Simulator of Polarized Hadronic INteract(X)ions. 

The features of Sphinx tt are similar to those of Sphinx, namely 

• Spin propagation. The spin information is propagated until the hard 
partonic interaction, i.e. the spin of the partons is taken into account in 
the initial state showering and hard partonic scattering. Final state in- 
teractions, such as final state showering and fragmentation, are treated 
as being spin- independent. 

• Polarized hard processes are described at leading order in the electro- 
magnetic and strong coupling constants. The processes currently im- 
plemented in the polarized mode are summarized in Table [^. Contact 
interactions are not available yet in the polarized mode. 



ISUB 


Process 


Comment 


11 


QiQj QiQj 


(anti-)quark - (anti-)quark scattering; 
annihilation is not included 


12 




annihilation nrocess 

CAjXXXXXXXX XCAj U X V-y XX K^-L K_j kj) 


13 


0:0; — >• 00 


annihilation into ffhion nair 

CAjXXXXXXXX XCAj U X V^XX X XX U X l_X V^XX Y-J C_AjXX 


14 


QiQi 91 


annihilation into gluon and prompt 7 


18 


QiQi 11 


annihilation into 7-pair 


28 


qi9 (li9 


(anti-) quark - gluon scattering 


29 


(li9 Qil 


prompt 7-production in (anti-)quark - 
gluon scattering 


53 


99 QkQk 


gluon fusion 


68 


99 99 


gluon - gluon scattering 


191 


QiQi fkfk 


annihilation into lepton-pair 
or quark - (anti-) quark pair 
(Drell-Yan process); this process is new 
and equivalent to the 7-piece of 
ISUB=1 in Pythia 



Table 1: List of processes implemented in the polarized mode 



• Polarized parton distributions. Two model parametrizations for 6q{x, Q"^) 
are provided, the inclusion of new ones is easy. 

• Documentation. Spin information has been added to the event listing. 

• Sphinx tt is equivalent to Pythia 5.6 in its unpolarized mode. 

Hereafter we do not try to explain the structure of Pythia 5.6 or the con- 
cepts that lead to the creation of Sphinx. Readers which are unfamiliar 
with these programs should consult @, 0. In the next section we will present 
the most important differences between Sphinx and Sphinx tt. Section ^ 
contains a detailed description of the modified subroutines as well as three 
tables with all new switches and internal variables. Finally an example main 
program is presented. 



2 From Sphinx to Sphinx TT: modifications 

The most notable difference between the scattering of two longitudinally po- 
larized nucleons and the scattering of two transversely polarized ones is that 
the symmetries in the cm frame are different. In the former case the system 
is invariant under rotations around the beam axis whereas in the latter case 
the transverse spin directions and the beam axis specify one or two planes in 



space. As a result, the rotational invariance is broken and the partonic cross 
sections now also may depend on the azimuthal scattering angle. 

In Sphinx tt the spins of the colliding nucleons are required to be either 
parallel or anti-parallel. Allowing for arbitrary transverse spins only changes 
the azimuthal dependence of the cross sections in a trivial way . 

Since the partonic cross sections in Pythia and Sphinx are independent 
of the azimuth it simply can be chosen according to a flat distribution. 
In comparison, the generation of the remaining kinematical variables is a 
complicated process. In order to make as few modifications as possible the 
event machinery in SPHINX TT was not rewritten but a trick was used to 
include the correct generation of the azimuth. 

During the initialization Pythia/ Sphinx searches for the absolute maxima 
of the various selected partonic subprocesses. If "by accident" the program 
finds a point in the kinematical space with a larger partonic cross section, 
the program adjusts the maximum value. For that reason, one can simply 
include the normal 0-dependent cross sections and all kinematical variables 
will be chosen correctly. 

This method only works for 2 —>■ 2-processes. Unfortunately, in Pythia / 
Sphinx the Drell-Yan process qiQi II was implemented as a 2 ^ 1-process 
with a successive decay of the gauge boson into a lepton pair. Furthermore 
the two steps were performed by different subroutines. It was impossible to 
implement transverse polarization here without violating the program struc- 
ture strongly and we simply included a new 2^2 subprocess for Drell-Yan 
events. For similar reasons the subprocesses which cover the QCD corrections 
to the simple Drell-Yan graph could not be made available. 

The modification of the initial state showering routine was relatively simple. 
Only the appropriate splitting functions had to be included and the handling 
of the Lorentz boosts required some changes in order to treat the azimuth 
correctly. 



3 Common Blocks and Subroutines 

In this section a detailed description of the modified subroutines is given. In 
order for this section to be self contained we present the differences between 
Pythia and Sphinx tt rather than giving a summary of the differences 
between Sphinx and Sphinx tt. The general structure of the code as well 
as the unchanged parameters and variables are not explained. However, it 
is shown in detail how the local polarization switch IPOL is implemented 
into the different subroutines. Only the polarized case (IPOL=l) will be 
discussed, because in the unpolarized case (IPOL=0) each subroutine works 
exactly as the corresponding Pythia 5.6 subroutine. The subroutines of 



Parameter 


Description 


Default 


MSTP(171) 


=0 
= 1 
=2 


beam polarization 
unpolarized 

polarization in +x direction 
polarization in —x direction 





MSTP(172) 


=0 
= 1 
=2 


target polarization 
unpolarized 

polarization in +x direction 
polarization in —x direction 





MSTP(175) 


use of polarized parton distributions in 
polarized initial state shower (MSTP(176)=l) 

=0: unpolarized distribution; for testing only; 
do not use! 

= 1: polarized 


1 


MSTP(177) 


=0 
= 1 
=2 

=3 


set of polarized parton distributions 5q{x, Q"^) 

used; in addition, one has to specify 

an unpolarized set as in standard Pythia 

6q{x, Q"^) = (no polarization) 

maximal polarization; 6q{x, Q"^) = q{x, Q"^) 

Model A: 5q{x,AGeV^) = Aq{x,AGeV^) 

+ correct DGLAP evolution 

set by Bartelski-Tatur is used 

same as before but Gehrmann-Stirling 

lo set C [H is used 





MSTP(178) 


percentage of beam and target polarization 





MSTP(179) 


switch off polarization temporarily in PYSIGH 

and PYSTFU resp. 
=0: no action 
=1: switch off polarization 




MSTP(180) 


mode selection (master switch) 
=0: unpolarized mode; this value overrides all 

other polarization switches 
= 1: polarized mode 





NSUB(ISUB) 


mode for subprocess ISUB 
=0: unpolarized treatment 
= 1: polarized treatment 





NSEL 


menu of polarized processes 
= 1: ISUB = 11,12,13,28,53,68 switched on 
= 10: ISUB = 14,18,29 switched on 






Table 2: Parameters controlling the polarized mode 



Variable 


Description 


Com. Block 


MINTOll) 


=1 
=2 


beam polarization 
unpolarized 

transverse polarization in +x direction 
transverse polarization in —x direction 


PYINTl 


MINT(312) 


=0 
=1 
=2 


target polarization 
unpolarized 

transverse polarization in +x direction 
transverse polarization in —x direction 


PYINTl 


MINT(313) 


=0 
=1 
=2 


polarization of shower initiator on beam side 
unpolarized 

transverse polarization in -\-x direction 
transverse polarization in —x direction 


PYINTl 


MINT(314) 


=0 
= 1 
=2 


polarization of shower initiator on target side 
unpolarized 

transverse polarization in +x direction 
transverse polarization in —x direction 


PYINTl 


MINT(315) 


=0 
= 1 
=2 


polarization of hard interacting parton 

on beam side 

unpolarized 

transverse polarization in +x direction 
transverse polarization in —x direction 


PYINTl 


MINT(316) 


=0 
= 1 
=2 


polarization of hard interacting parton 

on target side 

unpolarized 

transverse polarization in +x direction 
transverse polarization in —x direction 


PYINTl 


VINT (397) 


azimuthal angle of the hard scattering plane 
with respect to the x-axis before the initial state 
shower and primordial transverse momentum 
are generated; measured in cm-frame of hard 
interaction 


PYINTl 


VINT(398) 


same as VINT (397) after generation of initial 
state shower and prim, transverse momentum 


PYINTl 


VINT (399) 


scattering angle in the cm frame of the hard 

partonic interaction before the initial state 
shower and primordial transverse momentum 
are generated 


PYINTl 


VINT (400) 


same as VINT (399) after generation of initial 
state shower and prim, transverse momentum 


PYINTl 



Table 3: Internal variables storing polarization information 



Variable 



Description 



Com. Block 



ISIG(1000,6) 

ISIG(N,1) 

ISIG(N,2) 

ISIG(N,3) 

ISIG(N,4) 

ISIG(N,5) 

ISIG(N,6) 



hard scattering information of Nth channel 
particle code of Nth channel on beam side 
particle code of Nth channel on target side 
colour flow of Nth channel 
polarization of Nth channel on beam side 
polarization of Nth channel on target side 
not used 



PYINT3 



KD(I) 



polarization of Ith line in the 
event listing 
no polarization 

transverse polarization in +x direction 
transverse polarization in —x direction 



DPYPOL 



XSFX(2:-40:40,0:2) 



XSFX(JT,KFL,0) 
XSFX(JT,KFL,1) 

XSFX(JT,KFL,2) 



X times parton distribution for given 
X and of flavour KFL= -40:40 
and transvcrsity KFLD=0 : 2 on beam side 
(JT=1) and target side (JT=2) resp. 
unpolarized 

positive polarization with respect to 
transverse hadron polarization 
negative polarization with respect to 
transverse hadron polarization 



PYINT3 



Table 4: Internal variables storing polarization information 



Sphinx tt which we don't mention exphcitly remain unchanged with respect 
to Pythia 5.6. 

To incorporate polarization the following common blocks have either been 
enlarged, replacing the corresponding Pythia common blocks or have been 
added: 

• C0MM0N/PYINT3/XSFX(2,-40:40,0:2) ,ISIG(1000,6) , SIGH (1000) 

• C0MM0N/PYSUBS/MSEL,NSEL,MSUB(200) ,NSUB(200) , KFIN (2 , -40 : 40) , 
CKIN(200) 

• C0MM0N/DYP0L/KD(4000) 

Information about the new internal variables and enlarged arrays can be 
found in Table 0. 

MAIN PROGRAM 

Purpose: to set up the polarized event generation. The variables which 

have to be specified are listet in Table |^. 

Remarks: An example program is given in Section ^ 

SUBROUTINE PYINIT 
Purpose: 

• to display Sphinx tt header; 

• to check partially the availability of the desired polarization scenario, 
i.e. to check that the master switch for polarization MSTP(180) is 
set properly, that the selected partonic subprocesses can be treated 
polarized and to control and compose the polarization menue via NSEL; 

• to call DPLIST instead of LULIST (see below). 

New Parameters: MSTP(180), NSEL, NSUB(ISUB) 
Internal Polarization Switch: IPGL=MSTP(180) 

Remarks: If the chosen scenario is not allowed, the program is terminated 
with an appropriate error message. 

SUBROUTINE PYEVNT 
Purpose: 

• to start polarized event generation; 

• to call DPEDIT instead of LUEDIT (see below). 
New Parameters: MSTP(180) 

New Variables: VINT (397) , VINT(398) , VINT(399) , VINT(400) 
Internal Polarization Switch: IPGL=MSTP(180) 

Remarks: VINT (397) -VINT (400) can be used to measure the rotation of 



the hard scattering partonic system caused by the initial state shower and 
the primordial transverse momentum of the partons; if this effect is not small 
(check!) the transversity distribution is not used correctly in the hard scat- 
tering because the particle spins are no longer aligned along the x-axis; also, 
the implementation of the initial state shower is then no longer appropriate 
since at each partonic branching the azimuthal angles are chosen according to 
a flat distribution instead of the correct one. This latter shortcoming could 
be overcome if needed, but this would require major changes in the program. 

SUBROUTINE PYINKI (CHFRAM , CHBEAM , CHTARG , WIN) 
Purpose: 

• to check availability of the desired hadronic polarization scenario, i.e. 
to control that the selected hadron can be treated polarized; 

• to store the polarization of beam and target for the event listing in 
KD(1) and KD(2) and for internal use in MINT(311) and MINT(312). 

New Parameters: KD(I), MSTP(171), MSTP(172), MSTP(180) 
New Internal Variables: MINT(311), MINT(312) 
Internal Polarization Switch: IP0L=MSTP(180) 

Remarks: At present only nucleons and their antiparticles can be treated 
polarized. If a not allowed scenario has been chosen, the program is termi- 
nated with an appropriate error message. 

SUBROUTINE PYRAND 
Purpose: 

• to adapt PYRAND to the new enviroment - all relevant arrays which have 
been enlarged or added to the common blocks in other subroutines are 
modified here as well; 

• to extend event shape selection to incorporate transverse spin; 

• to store polarization of the partons entering the hard interaction in 
MINT(313) - MINT(316) (see Table |). 

New Parameters: MSTP(180), NSUB(ISUB) 

New internal variables: MINT(313) ,MINT(314) ,MINT(315) ,MINT(316) 
Internal polarization switch: IP0L=MSTP(180) xNSUB(ISUB) 
Remarks: Note that MINT(313)=MINT(315) and MINT(314)=MINT(316) 
but the values of MINT (313) and MINT (314) are changed later by the ini- 
tial state shower in PYSSPA. 

SUBROUTINE PYSCAT 
Purpose: 

• to adopt PYSCAT to the new environment (see PYRAND); 



• to store transverse spin information of the partons entering the hard 
interaction; 

• to add spin information to hnes 1, 2 and 5, 6 in the event hsting (see 
below) . 

New Parameters: KD(I), MSTP(180), NSUB(ISUB) 

New internal variables: MINT(315), MINT(316) 

Internal Polarization Switch: IP0L=MSTP(180) xNSUB(ISUB) 

SUBROUTINE PYSSPA ( IPUl , IPU2) 
Purpose: 

• to perform polarized initial state showering, transversity dependent 
DGLAP evolution equations are used in the backward evolution algo- 
rithm; 

• to enlarge all relevant arrays in an appropriate manner to incorporate 
polarization; 

• to check proper selection of the polarized initial state shower scenario, 
i.e. to control that MSTP(175) and MSTP(176) are set correctly; 

• to store the polarization of the initial state shower initiators (MINT (313) , 
MINT(314)). 

New Parameters: KD(I), MSTP(171), MSTP(172) , MSTP(175), 

MSTP(176), MSTPdSO), NSUB(ISUB) 

New Internal Variables: MINT (313) , MINT (314) , MINT (315) , MINT (316) 
Internal Polarization Switch: IP0L=MSTP(180) xNSUB(ISUB) 
Remarks: At the present stage only QCD shower can be treated polarized, 
QED showering has to be done in the unpolarized manner. The combination 
MSTP(175)=0 and MSTP(176)=1 allows to simulate polarized showering with 
the use of unpolarized parton distributions. This option is just for testing 
and should not be selected by the user! If MSTP(175) or MSTP(176) are set 
improperly, the program stops with an appropriate error message. The in- 
ternal variables MINT (313) and MINT (314) are changed to their final values 
in this subroutine. 

The handling of the Lorentz boosts in PYSSPA has been modified so that the 
azimuthal angle is propagated through the event. Nevertheless, at each par- 
tonic branching the azimuth of the decay plane with respect to the transverse 
spin of the decaying parton is distributed according to a flat distribution in- 
stead of the correct one. Only if the initial state shower does not cause a 
substantial rotation of the hard scattering partonic system this approxima- 
tion is justified. This shortcoming can be overcome if needed but would 
require substantial changes in the program. 

SUBROUTINE PYMULT(MMUL) 

Purpose: to switch off polarization in PYSIGH (set MSTP(179)=1 temporar- 
ily) when called from PYMULT. This is done because multiple interactions 



cannot be treated polarized at the moment. 
New Parameters: MSTP(179) 

SUBROUTINE PYREMN(IPU1 , IPU2) 
Purpose: 

• to adopt PYREMN to the new environment (see PYRAND) ; 

• to add spin information to hnes 3 and 4 in the event hsting (see below) . 
New Parameters: KD(I) 

Remarks: The generation of the primordial transverse momentum has been 
modified so that the azimuthal angles are now treated correctly. 

SUBROUTINE PYSIGH(NCHN,SIGS, IOC) 
Purpose: 

• to evaluate the transversity dependent hadronic cross sections by convo- 
lution of the spin dependent parton distributions with the transversity 
dependent partonic cross sections; 

• to supply the subroutine with the transversity dependent partonic cross 
sections. 

New Parameters: I0C,MSTP(171) ,MSTP(172) ,MSTP(179) ,MSTP(180) , 

NSUB(ISUB) 

Internal Polarization Switch: 
IP0L=MSTP(180) xNSUB(ISUB) x (1-MSTP(179) ) 
New Partonic Subprocess: ISUB=191 (g, + — > fkfk) 
Remarks: PYSIGH will always run in the unpolarizcd mode when it is 
called by PYMULT which sets temporarily MSTP(179)=1 in PYSIGH. When one 
evaluates the spin dependent hadronic cross sections one has to be aware 
that the hadrons are labelled according to their absolute spin whereas the 
partons are specified according to their relative spin with respect to the 
spin of the hadron. The parton distributions arc passed from PYSTFU to 
PYSIGH through the array XPQ (KFL , KFLD) (see below) and stored in the array 
XSFX(N,KFL,KFLD), where KFLD denotes the transversity. ISIG(N,I) con- 
tains the information about the Nth reaction channel of the chosen partonic 
subprocess. The new entries ISIG(N,4) and ISIG(N,5) specify the trans- 
verse spin of the partons at the beam and target side respectively. ISIG(N, 6) 
is reserved but not used at the moment. 

The cross sections in Pythia and Sphinx are independent of the azimuthal 
scattering angle which therefore can be easily chosen according to a flat dis- 
tribution. Contrarily, in SPHINX TT the cross section depends on the azimuth 
(f). Correct generation can be achieved by simply using the 0- dependent for- 
mulae of the partonic cross sections. Furthermore, during the initial search 
for the maximum of the partonic cross sections the azimuth must be fixed 
(I0C=0). During the normal event generation I0C=1 is set by PYRAND. No 



other changes are required to ensure correct generation of all kinematical 
variables. Contact interactions are not yet available for the polarized mode. 
The new partonic subprocess ISUB=191 that corresponds to ISUB=1 had to 
be included in PYSIGH in order to preserve the program structure. The im- 
portant difference is that no Z°-propagator is included in the new process 
and that ISUB=191 is a 2 — > 2-process rather than a2 1 one. 
One should be careful when using the subprocess ISUB=11 because here only 
interference graphs carry an asymmetry. These are only taken into account 
if one sets MSTP(34)=1. 

SUBROUTINE PYSTFU(KF,X,Q2,XPQ) 

Purpose: to evaluate the transversity dependent parton distributions for 
given flavour (KF), x (X), and (Q2) according to the selected parametriza- 
tions and models 

New Parameters: MSTP(177), MSTP(178), MSTP(179) , MSTP(180) , 
MSUB(ISUB) 

Internal Polarization Switch: IP0L=MSTP(180) x (1-MSTP(179) ) 
Remarks: Call of PYSTFU returns x times the parton distribution functions 
for given flavour, x, and for both spin orientations and an averaged (un- 
polarized) value. The values are stored in the array XPQ(KFL,KFLD) which 
has been enlarged from XPQ(-25:25) to XPQ(-25:25,0:2). XPQ(KFL,0) 
contains the unpolarized distributions, XPQ(KFL,1) distributions for posi- 
tive relative transversity and XPQ(KFL,2) distributions for negative relative 
transversity. The parton distributions are selected by switches described 
earlier in Table |^ The polarized distributions Q"^) are constructed from 

the unpolarized ones q{x, Q^) and the transversity distributions Sq{x, Q^) se- 
lected by MSTP(177) according to q-^i = |(g±(5g). The interface to the CERN 
parton distribution library has been updated so that Pdflib 4.0 and higher 
versions can be used with Sphinx tt. Two subroutines have been added 
to calculate the polarized distributions, namely BARTAT and GEHRMA. They 
require the data files bartat . dat and gehrmann . dat which are supplied with 
the program and must be visible to the FORTRAN open statement. Therefore 
the paths of these files have to be modified appropriatly in the subroutines 
mentioned above. When q and Sq are combined to compute the distributions 
for fixed transversity an unitarity check is performed — if one of the result- 
ing values is negative this value is set to zero and the value for the other 
spin direction is set equal to the unpolarized one. Only polarized proton 
parametrizations are implemented. Neutron distributions are obtained by 
SU(2) symmetry. Charge conjugation is used to describe the corresponding 
antiparticles. 

SUBROUTINE PYSTPR(X,Q2,XPPR) 

Purpose: to calculate the unpolarized parton distributions of the proton; 
to update the implementation of the Cteq2 parton distribution functions 
Remark: the Interface to PYCTQ2 was taken from Pythia 5.7 



FUNCTION PYCTQ2 (Iset, Iprt , X, Q) 

Purpose: to give the revised Cteq2 parton distribution sets with extended 
range in parametrized form. 

Remark: This function was taken from Pythia 5.7. 
SUBROUTINE BARTAT (X , Q2 , UPV , DNV , EM , DELTA) 

Purpose: to return x times the polarized parton distributions evaluated at 
given X (X) and (Q2) for model 1. UPV, DNV, EM, DELTA contain 6uy, 
6d^, 6M and 66 ^ repectively. 

Remarks: We use 6q{x, 4 GeV^) = Aq{x, 4 GeV^) for lack of any knowledge 
of 6q{x,Q'^). One actually has to expect that 6q and Aq differ substan- 
tially because their evolution equations are different. For Ag(a;, 4GeV^) the 
Bartelski-Tatur parametrization is used. The (5g-values for higher 
were obtained by correct DGLAP evolution. BARTAT requires the data file 
bartat .dat. 

SUBROUTINE GEHRMA (X , Q2 , UPV , DNV , EM) 

Purpose: to return x times the polarized parton distributions evaluated at 
given X (X) and Q"^ (Q2) for model 2. UPV, DNV, EM contain 6uy, 6dy and 6s 
repectively. The sea is assumed to be SU(3) symmetric. 
Remarks: In this model 6q{x, AGeV'^) = Ag(x,4GeV^) is assumed for lack 
of any knowledge of 6q. The Germann-Stirling Set C is used for Aq. The 
5g-values for higher were obtained by correct DGLAP evolution. GEHRMA 
requires the data file gehrmann.dat. 

SUBROUT I NE DPL I ST ( ML I ST ) 

Purpose: to display the polarizations of the particle in the event listing. 
New Parameters: KD(I) 

Remarks: DPLIST is a modification of the Jetset subroutine LULIST. It is 
changed to display the polarization in the final listing. The sign displayed 
just behind the particle code denotes polarization with respect to the 
When the sign is missing the particle has been treated as being unpolarized. 
The information is taken from the vector KD(I) and printed using the symbols 

= ' ' : no polarization (KD(I)=0) 

= 'A' : transverse polarization in +x direction (KD(I)=l) 
= 'v' : transverse polarization in —x direction (KD(I)=2) 

SUBROUTINE DPEDIT(MEDIT) 

Purpose: to properly compress the vector KD(I) which contains the polar- 
ization information. 
New Parameters: KD(I) 

Remarks: DPEDIT is a modification of the Jetset subroutine LUEDIT. 



SUBROUTINE PDFSET (FARM, VALUE) 

Purpose: new dummy routine for new Pdflib interface; has to be removed 
if Pdflib is used. 

SUBROUTINE STRUCTM (XX , QQ , UPV , DNV , USEA , DSEA , STR , CHM , BOT , TOP , GLU) 
Purpose: new dummy routine for new Pdflib interface; has to be removed 
if Pdflib is used. 



4 The Main Program 



4.1 The code 



PROGRAM EXAMPLE 

C Example of a Main Program for event generation 

C in proton-proton scattering with transverse polarization 

C 

C One beam is polarized in +x direction, the other one in 

C -X direction 

C 

C This program has to be linked with 

C the programs SPHINX TT and JETSET 7.3, 

C the data files BARTAT.DAT and GEHSTI.DAT, and 

C the CERN Libraries 

C COMMON BLOCKS of SPHINX TT for event generation 

C0MM0N/LUDAT1/MSTU(200) ,PARU(200) ,MSTJ(200) ,PARJ(200) 
C0MM0N/LUJETS/N,K(4000,5) ,P(4000,5) ,V(4000,5) 
C0MM0N/PYSUBS/MSEL,NSEL,MSUB(200) ,NSUB(200) , 
& KFIN (2 , -40 : 40) , CKIN (200) 

C0MM0N/PYPARS/MSTP(200) ,PARP(200) ,MSTI(200) , PARI (200) 
C0MM0N/PYINT5/NGEN (0 : 200 , 3) , XSEC (0 : 200 , 3) 
COMMON/PAWC/HBOOK (10000) 

DIMENSION R0B0(3) 

C============================================================= 

C polarization set up 

C============================================================= 



C polarized simulation 



MSTP(180)=1 



C beam polarized in +x direction 
MSTP(171)=1 

C "target" polarized in -x direction 
MSTP(172)=2 

C "fully" polarized parton distributions 
MSTP(177)=1 

C the degree of polarization of each beam is only 70% 
MSTP(178)=70 

C polarized initial state shower 
MSTP(176)=1 

C============================================================= 

C event set-up 

C============================================================= 

C select hard subprocesses a la carte 
MSEL=0 

C choose "Drell-Yan" q q~ — > gamma — > f f~ 
MSUB(191)=1 

C switch on polarization for 191 
NSUB(191)=1 

C switch off multiple interactions 
MSTP(81)=0 
MSTP(131)=0 

C number of events to be generated 
NEVENT=1000 

C============================================================= 

C start 

C============================================================= 

C initialize 

CALL PYINITCCMS' , 'p' , 'p' ,200.) 

C book histograms 

CALL HLIMIT (10000) 

CALL HBQOKKl, 'azimuthal distrib. ' ,20,-3. 1416,3. 1416,0. ) 



C event loop 

DO 200 I=1,NEVENT 



CALL PYEVNT 

C list first few events 

IF(I.LE.6) CALL DPLIST(l) 

C detect myons and electrons 

IF(ABS(K(7,2)) .EQ.13.0R.ABS(K(7,2)) .EQ.ll) THEN 

C boost to di-lepton rest-frame 

DO 100 J=l,3 

R0B0(J)=(P(7,J)+P(8,J))/(P(7,4)+P(8,4)) 
100 CONTINUE 

ROBOT=ROBO (1) **2+R0B0 (2) **2+R0B0 (3) **2 
IF (ROBOT. GE. 0.999999) THEN 

R0B0T=1 . 0000 1*SQRT (ROBOT) 

R0B0(1)=R0B0(1) /ROBOT 

ROBO (2) =ROBO (2) /ROBOT 

ROBO (3) =ROBO (3) /ROBOT 
ENDIF 

CALL LUDBRB(1,N,0. ,0. , -DBLE(ROBO (1) ) , 
& -DBLE(R0B0(2)) ,-DBLE(R0B0(3))) 

C calculate azimuthal angle 
PHI=PLU(7,15) 

C fill histogram 

CALL HFILL(1,PHI,0. ,1.) 
ENDIF 

200 CONTINUE 

C print cross section and statistics 
CALL PYSTAT(l) 

CALL HRPUT (0 , ' example . hbk ' , ' n ' ) 



END 



4.2 The Event Listing 



SPHINX TT version 1.00 
** Last date of change: 10 Oct 1996 ** 



The Lund Monte Carlo - PYTHIA version 5.6 
** Last date of change: 3 Apr 1992 ** 



The Lund Monte Carlo - JETSET version 7.3 
** Last date of change: 10 Mar 1992 ** 

1************* PYINIT: initialization of PYTHIA routines ******** 



I I 

I PYTHIA will be initialized for a p on p collider I 

I at 200.000 GeV center-of-mass energy I 

I I 



** PYMAXI: summary of differential cross-section maximum search * 



I II 
I ISUB Subprocess name I Maximum value I 

I II 



I II 
I 191 f + f~ -> f + f~ I 6.5589E-05 I 

I II 



*************** PYINIT: initialization completed **************** 



Event listing (summary) 



I particle/ jet KF orig p_x 



p_y p_z E m 



1 


!p+! 


2212-^ 





0, 


.000 


0, 


.000 


99, 


.996 


100, 


,000 


0, 


,938 


2 


!p+! 


2212V 





0, 


,000 


0, 


,000 


-99, 


.995 


100. 


,000 


0, 


,938 


3 


!u! 


2" 


1 


-0, 


,414 


-0, 


.069 


25, 


.831 


25, 


,835 


0, 


,000 


4 


!u~! 


-2v 


2 


-0. 


,207 


-0, 


.420 


0, 


.321 


0, 


.568 


0, 


,000 


5 


!u! 


2" 


3 


-0. 


,414 


-0, 


.069 


25, 


.831 


25, 


.835 


0, 


,000 


6 


!u~! 


-2v 


4 


-0, 


,207 


-0, 


,420 


0, 


.321 


0, 


.568 


0, 


.000 


7 


!u! 


2 





-1. 


,025 


0, 


,598 


23, 


. 197 


23, 


.227 


0, 


,006 


8 


!u~! 


-2 





0, 


,404 


-1, 


,087 


2, 


.956 


3, 


.175 


0, 


,006 



Warning: maximum violated by 1.017E+00 in event 10 
ISUB = 191; Point of violation: 

tau=4.112E-04, y*=-l .309E+00, cthe=-0. 6715451, tau'=0.000E+00, 
phi = 1.299 

XSEC(191,1) increased to 6.673E-05 

Warning: maximum violated by l.OllE+00 in event 82 
ISUB = 191; Point of violation: 

tau =3.328E-04, y*=-8.523E-01, cthe=0. 8431695, tau'=0.000E+00, 
phi = 1.961 

XSEC(191,1) increased to 6.745E-05 

Warning: maximum violated by 1.039E+00 in event 566 
ISUB = 191; Point of violation: 

tau =3.455E-04, y*=-9 . 051E-01 , cthe=0 . 8579352, tau'=0.000E+00, 
phi = 1.735 

XSEC(191,1) increased to 7.010E-05 

1*** PYSTAT: Statistics on Number of Events and Cross-sections *** 



I III 

I Subprocess I Number of points I Sigma I 

I III 

I 1 1 (mb) I 

I III 

I N:o Type I Generated Tried I I 

I III 



I III 

I All included subprocesses I 1000 4831 I 1.481E-05 I 

I 191 f + f~ -> f + f~ I 1000 4831 I 1.481E-05 I 

I III 



*** Fraction of events that fail fragmentation cuts = 0.00000 *** 
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